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Abstract. We address consequences of strong tensor and weak spin-orbit terms in the local energy density 
functional, resulting from fits to the /5/2 — /7/2 splittings in '^''Ca, ''^Ca, and ^^Ni. In this study, we focus 
on nuclear binding energies. In particular, we show that the tensor contribution to the binding energies 
exhibits interesting topological features closely resembling that of the shell-correction. We demonstrate that 
in the extreme single-particle scenario at spherical shape, the tensor contribution shows tensorial magic 
numbers equal to A(Z)=14, 32, 56, and 90, and that this structure is smeared out due to configuration 
mixing caused by pairing correlations and migration of proton/neutron sub-shells with neutron/proton 
shell filling. Based on a specific Skyrme-type functional SLy4T, we show that the proton tensorial magic 
numbers shift with increasing neutron excess to Z=14, 28, and 50. 

PACS. 21.60.Jz Nuclear Density Functional Theory and extensions (includes Hartree-Fock and random- 
phase approximations) 21.60.-n Nuclear structure models and methods 



1 Introduction 

Density functional theory (DFT) is a universal ab ini- 
tio approach designed and used to calculate properties 
of electronic systems entrapped in the external Coulomb 
field of nuclei. It has been successfully applied to atoms, 
molecules, or condensed media. Universality of the DFT 
means independence of a functional form of the shape of 
external one-body potential holding the electronic system 
together. The existence of such a universal and, in prin- 
ciple, exact density functional describing ground-states of 
externally bound fermionic systems is guaranteed by the 
Hohenberg-Kohn [1] and Kohn-Sham [2] (HKS) theorems. 

Generalization of the DFT theory to self-bound sys- 
tems like atomic nuclei encounters problems associated 
with intrinsic rather than laboratory density which char- 
acterizes the atomic nuclei, see Ref. [3l|4j[5]. In spite of 
that, the HKS theorem have strongly influenced the way 
of thinking in the field of nuclear structure. Nowadays, 
the nuclear structure theorists employ the functionals that 
are treated as separate entities, independently, to a large 
extent, of the underlying effective nucleon-nucleon (NN) 
interactions like, for example, the local Skyrme interac- 
tion [6]. Free parameters of these functionals are directly 
adjusted to fit empirical data. There are also attempts 
to enrich their functional form as compared to the form 
resulting from the mean-field (MF) averaging of the ef- 
fective Skyrme interaction [7j, which are motivated by a 
rather mediocre performance of the conventional Skyrme- 
type functionals, see for example Refs. fS'j'S]. 



Adequateness of the fitting strategy - that is, the choice 
of the data set - is a key factor determining performance 
of the nuclear energy density functional (EDF) method. 
In this work we explore the Skyrme-type local EDF ap- 
proach to nuclear structure and focus on the spin-orbit 
(SO) and tensor parts of the functional. Throughout the 
years, not much attention was paid to, in particular, the 
tensor part, mainly due to the lack of clear experimental 
data constraining the strength of this part of the EDF. 
Hence, the tensor term in the existing Skyrme function- 
als is either trivially set to zero by hand, see the review 
in Ref. [TD], or is a result of a global fit to bulk nuclear 
properties [TT] . 

Recent revival of interest in the tensor term was trig- 
gered by systematic observation of non-conventional shell 
evolution in isotopic chains of light nuclei far from stabil- 
ity, including new magic shell-gap opening at iV=32, see 
for example Refs. [T21IT3] . Otsuka and collaborators associ- 
ated these effects with the two-body tensor interaction [HI 
115] . This interpretation was soon confirmed within the lo- 
cal Skyrme-type EDF models [l6l[T7l[T8l[T9l[20l[2T] . Inclu- 
sion of single-particle (s.p.) energies in the fitting data sets 
appears to lead to the tensor coupling constants 1 171! 201 1211 
[22 l [23 j [24] which are at variance with bulk fits, see Fig. [1] 
and extensive discussion in Ref. [11 . 

The aim of this work is to look into consequences of 
strong attractive isoscalar and isovector tensor fields and 
weak SO fields resulting from the fitting method proposed 
by our group [55]. In particular, we study such conse- 
quences for the nuclear binding energies. The paper is 
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Fig. 1. Proton Cp = {Ci - Ci^)/2 versus neutron Ci = (Cq + 
C()/2 tensor coupling constants resulting form fits to: bulk 
nuclear properties (black dots) and the s.p. levels and the SO 
splittings (triangles). Open triangles represent fits of Ref. [171 
I20II21| . Black triangles mark our results | ,22i.23..24) from the fit 
to the /7/2 — fs,/2 SO splittings. Shaded area represents the 
parameters established by Brink, Stancu and Flocard (BSF) 
in their seminal paper ^25^ 



organized as follows. In Sec. [21 theoretical framework is 
briefly outlined. In Sec. [3l the procedure used to fit the 
tensor and SO coupling constants is discussed. In Sec. [H 
numerical results showing tensor energy contribution to 
the total binding energy, followed by a discussion of ten- 
sorial magic structure, is presented. The paper is summa- 
rized in Sec. [5] 



2 Tensor and spin-orbit parts of the local 
nuclear energy density functional 

In this work we consider the local EDF H{r) of the Skyrme- 
type. It consists of a kinetic energy and a sum of isoscalar 
{t=0) and isovector {t=l) potential energy terms: 



where 



(1) 



(2) 



(3) 



with the density-dependent primary coupling constants 
Cf [po] and [po]- The potential energy terms are bihn- 
ear forms of either time-even (p, r, J) or time-odd (s, j. 



T) densities and their derivatives, see, e.g.. Ref. [10] for 
details. The density denotes the vector part of the spin- 
current tensor, 3t,\ = I]^^. <^x^L^It,^L^- 
In this work we focus on the tensor, 

= C^J§ + C/J2, (4) 

and the SO terms, 

T^so = Co^^poV- Jo + C7VV-Ji. (5) 

In the limit of spherical symmetry, the vector part Jt = 
Jt(r)er is the only non- vanishing part of the tensor density 
Jp^. Hence, in this limit, the tensor part of the functional 
dll reduces to: 



(6) 



By performing variation of the functional with respect to 
Jt{r) one obtains the one-body SO potential: 



1 

2^ 



J dpt 



C/Jt(r) L-S 



(7) 



which is composed of two terms. The first term is coming 
form the SO term in the functional, Eq. ([5]). It is propor- 
tional to the radial derivative of the particle density and is 
relatively slowly varying with N and Z. The second com- 
ponent is due to the tensor term Q. It is proportional 
to the SO density Jt(r) which is strongly shell- filling de- 
pendent. Indeed, in the spherical symmetry limit, the SO 
vector density can be written as [26] : 



J(r) 



njl 



j{j + 1) - l{l + 1) 



€Jr), (8) 



where v^ji and i'njii^) occupation probabilities and 
radial wave functions, respectively, of states with given 
quantum numbers. If both SO partners = I ±1/2 are 
fully occupied, i.e., when the system is spin-saturated (SS) 
the J(r) density vanishes. Examples of the SS systems in- 
clude ^^O, ""^Ca, or ^''Zr at spherical shape. Most of the 
nuclei are spin- unsaturated (SUS). The SO vector density 
reaches its maximum when one (or more) of the SO part- 
ners is fully occupied while the other one is completely 
empty. 



3 Fitting procedure 



Fitting procedure used to constrain the coupling constants 
C/ and C^"^ was described in detail in Ref. [21] and we 
only recall it here very briefly. The idea is to reproduce 
experimental /7/2 — /5/2 SO splittings in three key nuclei: 
^"Ca, ''^Ca, and ^^Ni. Since ^"Ca is, as discussed above, 
a SS system, the conventional SO term of Eq. ^ is the 
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Fig. 2. Empirical (horizontal lines) and theoretical (inclined 
lines) /7/2-/5/2 splittings in *°Ca, ***Ca, and ^^Ni as a function 
of C^-^ /C^\ShyA) (upper panel), Ci (middle panel), and C/ 
(lower panel). Solid an dashed lines represent neutron an pro- 
ton splittings, respectively. The theoretical results are obtained 
by modifying the SO and tensor strengths in the SLy4 func- 
tional. Empirical data are taken from [27) . See text for further 
details. 



only source of the SO splitting. Hence, this nucleus is used 
to set the isoscalar strength of the SO term Cq'"^. Having 
set Co''', we next constrain the Cg strength by using the 
h/2 — SO splitting in the isoscalar (A^ = Z) SUS 
nucleus ^^Ni. Finally, we move to ^^Ca, where protons 
and neutrons constitute a SS and SUS system, respec- 
tively. This nucleus is used to fit the isovector coupling 
constants or, more precisely, to fit C7"^, because the ra- 
tio Cq''^ /CY'^ is kept equal to the value characteristic for 
the given Skyrme parametrization. There is one piece of 
data on the /7/2 — /5/2 SO spHttings, preferably in ^^Ni 
or ^^Ni, which is badly needed to fit the tensor and SO 



dent on the initial parameterization. This is illustrated in 
Fig. [1] where different functionals modified according to 
our prescription, see Rcfs. [22, 23, 24J, are collected. They 
are labeled by a subscript T following the force acronym 
and marked by black triangles. 



4 Contribution from the tensor terms to the 
binding energy 

Fig. [3] shows the contribution to the total nuclear bind- 
ing energy due to the tensor term, calculated by using the 
spherical Hatree-Fock-Bogolubov (HFB) code HFBRAD 
[29] with the SLjAt functional. Contributions due to the 
isovector and isoscalar parts are depicted separately in the 
upper and middle panels, respectively. The total contribu- 
tion, Bt{N, Z), is shown in the lowest panel of the figure. 

From these results one can see that the isovector com- 
ponent is rather weak. Hence, the topology of the total 
contribution to the energy is mostly determined by the 
isoscalar term that shows a strong shell dependence. Fol- 
lowing the argumentation presented in Sec.[2l the strongest 
tensor effects are expected to appear for N{Z) = 14, 32, 
56, and 90. They correspond to nucleons filling the 1^5/2 1 
1/7/2 ® 2p3/2, I59/2 © 2^5/2 and I/111/2 ® 2/7/2 shells, re- 
spectively, which creates a maximum SUS filling. Hence, 
these numbers can be regarded as tensorial magic num- 
bers. 

However, as it is seen from Fig.dS]), the Bt{N, Z) does 
not follow the expected pattern exactly. This is due to 
(i) pairing-induced configuration mixing and {iij changes 
in the s.p. ordering of levels caused by the combination 
of strong attractive tensor fields and strongly reduced SO 
field. Two such situations are visible in Fig.[3l For iV < 30 
the tensor contribution is, as expected, largest for 32 Ge. 
For 40 < < 50, however, the minimum on the plot is 
shifted toward the 28^1 isotopes, which suggests a change 
in the order of the IP3/2 and I/5/2 proton sub-shells with 
increasing neutron excess. The figure also indicates that 
on the proton side, Z = 50 rather than Z = 56 is the ten- 
sorial magic number. Again, this suggests that the I57/2 
proton sub-shell is filled before 1^5/2- Consequently, the 
tensorial magic numbers may slightly differ for neutrons 
{N = f 4, 32, and 56) and for protons {Z = f 4, 28, and 50). 
This effect, however, may strongly depend upon a rather 
delicate balance between the tensor and SO strengths and 
needs to be studied in detail. 



terms uniquely. 

The procedure outlined above is illustrated in Fig. [2] 
for the case of the SLy4 functional [28 but it is quali- 5 Summary 
tatively independent of the initial parameterization. As 
shown in the figure, a good agreement with empirical data 
requires, for this low-effective-mass force, circa 35% reduc- 
tion of C^^ as compared to the original value, a large at- 
tractive isoscalar tensor coupling constant C^'^ of about 
-45MeVfm^ and CY'' of about -TOMeVfm^ ft appears 
that the resulting tensor coupling constants C^'^ (as well 
as the SO strengths C^) are, to large extent, indepen- 



A new strategy of fitting the coupling constants of the nu- 
clear energy density functional was recently proposed by 
our group [52] . It is based on a fit of the isoscalar spin-orbit 
and both isoscalar and isovector tensor coupling constants 
directly to the the - /7/2 SO splittings in ^^Ca,, ^^Ni, 
and ''^Ca. Our results show that drastic changes in the 
isoscalar SO strength and the tensor coupling constants 
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Fig. 3. The isovector (upper), the isoscalar (middle) and the 
total (lower) tensor contribution to the nuclear binding en- 
ergy. The calculations were done using the SLy4T interaction 
in the particle-hole channel and the volume-(5 interaction in 
the particle-particle channel with spherical symmetry assumed. 
Vertical and horizontal lines indicate the single-particle tenso- 
rial magic numbers at spherical shape. See text for further 
details. 



are required as compared to the commonly accepted val- 
ues 

This work briefly addresses the consequences of strong 
attractive tensor and weak SO fields on binding energies. 
In particular, a contribution to the binding energy due 
to the tensor interaction is calculated. It shows a generic 
pattern closely resembling that of the shell-correction. The 
tensorial magic numbers are shifted up relatively to the 
standard magic numbers, towards N{Z)=\A, 32, 56, or 90, 
which, in the extreme s.p. scenario at the spherical shape, 
correspond to the maximum spin-orbit asymmetry in the 
1^^5/2, 1/7/2 ©2p3/2, lff9/2©2d5/2, and I/111/2® 2/7/2 con- 
figurations, respectively. It is shown that these numbers 
are smeared out by the pairing effects and shifted in the 



case of protonic tensorial magic numbers by changes in the 
sub-shell ordering. It is also shown that strong attractive 
tensor interaction may give rise to an increased stability 
of nuclear binding at the drip lines, in particular around 
Z « 14, TV « 32, Z « 28, N « 56, and Z « 50, TV « 90. 

This work was supported in part by the Polish Ministry 
of Science under Contract No. N N202 328234, by the 
Academy of Finland and University of Jyvaskyla within 
the FIDIPRO programme, and by the Swedish Research 
Council. 
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